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Le but de ce cours est d ' apprendre a programmer en fortran 90. II est prevu 
pour 8 seances d'lh20, accompagne de 8 seances de quatre heures de TD-TP 
sur machine. 

J'ai essaye de le presenter non pas comme une collection exhaustive de com- 
mandes (ce qui permet difficilement de pouvoir programmer avant d'avoir fini 
le cours) mais plutot comme une progression en trois phases. La premiere par- 
tie devrait permettre de pouvoir ecrire vos premiers programmes apres une ou 
deux seances seulement. La deuxieme introduit certains concepts plus complexes 
(sous-programmes, types derivees, allocation dynamique). Les troisieme et qua- 
trieme parties mettent l'accent sur les apports reellement nouveaux de fortran 
90 par rapport au vieux fortran 77 (modules, interfaces generiques, surcharge 
d'operateurs, pointeurs). 

Nombre de commandes presentees ici sont expliquees avec beaucoup plus de 
precision dans divers cours et manuels de reference donnes dans la bibliographic 
Ces references sont en acces libre et gratuit sur internet. N'hesitez pas a les 
consulter des que vous en avez besoin. 

Le langage fortran evolue regulierement , mais la veritable evolution a ete faite 
en passant du fortran 77 au fortran 90. Les nouvelles normes fortran 95 (puis 
fortran 2003, fortran 2008) n'apportent que des fonctionnalites nouvelles. Par 
consequent, tout ce qui est dit ici devrait rester valable avec les normes futures. 

Enfin, je veux insister sur le fait que quel que soit le cours que vous utiliserez, il 
n'y en fait qu'un seul moyen d'apprendre a programmer : c'est de programmer 
soi-meme (de la meme facon que l'on n'apprend bien une langue etrangere qu'en 
la parlant soi-meme). Done usez et abusez de votre temps en salle machine (le 
temps des seances de TP n'y suffira pas), c'est la garantie du succes. 
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1 Introduction 

1.1 Qu'est-ce qu'un programme infor- 
matique ? 

tout commence par un algorithme : suite d 'operations 
donnant un resultat final 
exemple : 

approximation de \/2 par methode de Newton 
uo donne 

calculer Un+i = — H — — , pour n = a 10. 

2 Un 

- ces operations sont ecrites dans un fichier, avec un 
langage informatique (pour nous, le fortran 90 (f90)) 
—> ce fichier constitue le programme ; 
— ► c'est une suite d'operations sur des variables, et 

event uellement des echanges de donnees entre uti- 

lisateur et ordinateur ; 
—> les variables sont des noms correspondant a de 

la place memoire dans l'ordinateur, utilisee pour 

stocker les donnees et les manipuler. 
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> exemple : approximation de y/2 



program racine 




implicit none 




/ variables 

integer n 
real u 




/ initialisation 

u = 1.0 




/ b o u c I e 

do n = l,10 

u=u/2.0 + 1.0/u 
end do 




/ a f f i ch a g e 

pri n t * , ' approx . de sqrt (2) : 




end program racine 
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ce programme n'est pas utilisable tel quel : il faut 
le rendre comprehensible par l'ordinateur, c.-a-d. le 
traduire : le logiciel traducteur est appele compila- 
teur, la traduction est appelee compilation. 
exemple : 



gfortran -o prog prog.f90 



prog. f 90 : fichier contenant le programme ecrit en 



f90 



prog I : traduction (appele fichier executable) 



le fichier executable permet d'executer les instruc- 
tion programmees. 
exemple : 

./prog 
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1.2 Le langage F90 : generalities 

- ensemble de regies syntaxiques et grammaticales et 
de mots permettant d'ecrire des operations mathe- 
matiques, de communiquer des informations entre 
l'ordinateur (memoire), les fichiers (disque dur) et 
l'utilisateur (clavier/ecran). 

l'essentiel du cours va porter sur l'apprentissage de 
ce langage : 

— ► en commencant par une vue rapide et simplifiee 
permettant de pouvoir programmer apres une ou 
deux seances ; 

— ► puis en decouvrant ensuite les aspects plus com- 
plexes du langage. 

- conventions typographiques utilisees dans ce docu- 
ment : en gras pour les mots-clefs, entre crochets 

pour les instructions optionnelles, en italique pour 
les commandes unix, en teletype pour les lignes 
fortran. 
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- historique : premier compilateur (langage de haut 
niveau) en 1957, evolutions constantes depuis (voir 
wikipedia) : 

- 1978 : Fortran 77, tres repandu aujourd'hui encore dans les 
gros codes de calculs 

- 1990 : Fortran 90, fortran devient un langage moderne (format 
libre, programmation modulaire) 

- fortran 95, 2000 (programmation orientee objet), 2004, ... 

F90 ou C++ pour le calcul scientifique : la majeure 
partie des gros codes de calcul industriel sont ecrits 
en fortran ... 

- Attention aux effets de mode ("Fortran est depasse", 
"vous programmez encore en fortran ? v , "il faut pro- 
grammer en C-h-h v ) : 

- fortran est bien plus simple d'utilisation 

- les concepts de C++ ne sont guere utiles pour le calcul scien- 
tifique en general 

- les aspects programmation "orientee objet" du C++ sont main- 
tenant utilisables en fortran (depuis fortran 90 et fortran 2003) 

- de nouveaux compilateurs fortran sont developpes pour les 
nouvelles architectures (fortran pour les GPU par exemple). 

- attention aux commentaires de ceux qui ne connaissent que le 
C (et de ceux qui ne connaissent que le fortran!) 

- voir le blog d'un physicien qui connait bien les deux langages : 
http : //qd . typepad . com/ 4/ 2005/ 08/ c_vs_f ortran . html 
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2 Premiers pas en F90 

2.1 Exemple de programme 

2.2 structure 

program nom_du_programme 

declaration des variables 
instructions 

end program nom_du_programme 



2.3 Les variables 

- variables : noms (chaines de caracteres alphanumeriques 
plus quelques autres, cf. section 2.4) permettant de 
manipuler des donnees en memoire. 

- operation essentielle : V affectation (symbole =) 
exemple : variable n (entier) : 

n=2 / — p r end la v a I eur 2 
n=n+ 1 / — au g m en t e de 1 

ce n'est pas une egalite mais une operation : la valeur 
de n stockee en memoire est remplacee par elle-meme 
plus 1. 
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- le bloc declaration des variables sert a indiquer a la 
machine quel est le type de chaque variable utilisee 
afin de reserver la place en memoire necessaire. 
exemple : 



declaration 


signification 


ex. d'affectation 


integer : : n, i 


entier 


n=10 


real : : x,y 


reel 


x=1.0e-5 






y=1.0 


complex : : z 


complexe 


z=cmplx (1.0,2. 5e3) 


character : : c 


caractere 


c= ; d ; 


logical : : b 


booleen 


b= . true . 



- declaration implicite : bien que f90 ne l'exige pas, 
il est prudent de declarer toutes les variables. Pour 
eviter les oublis, placer l'instruction implicit none 
au debut du bloc declaration. Les oublis sont alors 
detectes a la compilation. 
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- les tableaux : tous les types precedents peuvent etre 
structures en tableaux avec l'attribut dimension 




integer, dimension ( —2:2) dim 
real, dimension ( 1 : 3 , : 4 , 1 : 2 ) tab 
character, dimension ( 1 : 3 , 4 : 6 ) ten 




—> chaque argument de dimension indique les bornes 
de la dimension correspondante par un intervalle en- 
tier mimmax 

—> acces aux composantes : 




dim(0) = 10 

tab (2 ,3 ,20) = 2.3e-5 
tch(l ,5)= 'a ' 
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2.4 les identificateurs 

- caracteres du clavier permettant de nommer une va- 
riable ou un autre element du programme ; 

- suite de 1 a 31 caractere alphanumeriques parmi : 
lettres sans accents majuscules et minuscules, chiffres, 
caractere _ (souligne) ; 

- le l er caractere doit etre une lettre ; 

- les majuscules et minuscules ne sont pas different iees ; 

- exemple : 

—> les identificateurs suivants sont valides : 

const ant e_gaz 

pi2 

Rk54 

— ► les identificateurs suivants ne sont pas valides : 

accentue 
avec espace 

Il_y_a_plus_de_trente_et_un_caracteres 
_souligne_devant 
1 _chif f re.devant 
nom#alphanumer ique 
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2.5 Continuation de ligne : caractere & 

- une ligne d'instructions comporte au maximum 132 
caracteres 

- si plus de 132 =4> message d'erreur : 

write (*,*) 'Au clair de la lune, mon ami Pierrot. Prete moi ta plume, pour ecrire un mot. 
Ma chandelle est morte, je n ai plus de feu. Ouvre moi la porte pour 1 amour de dieu' 
1 

Error : Unterminated character constant beginning at (1) 



- solution : couper la ligne en morceaux. Placer & en fin 
de premiere ligne et & au debut de la ligne suivante : 



write ( * , * ) 


' Au clair de la lune , mon ami 


Pierrot . & 




& 


Prete moi ta plume , pour e 


c r i r e un mot 


. & 


& 


Ma chandelle est morte, je n 


ai plus de 


& 


& 


feu . Ouvre moi la porte pour 


1 amour de d 


i e u ' 
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2.6 Les comment aires : symbole ! 

- tout ce qui est a droite du symbole ! est ignore par 
le compilateur 

- cela permet d'ecrire des commentaires dans le pro- 
gramme, c.-a-d. du texte qui en explique le fonction- 
nement 

- le ! peut etre place n'importe ou dans la ligne (au 
debut, au milieu, a la fin), et a n'importe quel endroit 
du programme 

- exemple : voir le programme page 30. 
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2.7 Les structures de controle 

Ce sont des instructions qui permettent de faire des 
iterations ou des tests sur des variables. 

2.7.1 Le test if 

- exemple : calcul de la valeur y = { xl ^ x S1X ^° 

^ & v sinon 



if ( x / = 0.0) then 

y = x*log(x) 
else 

y = o.o 

end i f 



syntaxe generale : 



nom : ] if (expression logique) then 

instructions 
elseif (expression logique) then 

instructions 

else 

instructions 
end if [nom 



ou nom permet de donner un nom a la structure pour 
une meilleure lisibilite du programme (par exemple 
en cas de if imbriques) 
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- expression logique est une combinaison d'un ou 
plusieurs tests sur des variables, ainsi que des trois 
operateurs logiques 



test 


signification 




egal 




dirterent 


> 


superieur strict 


< 


inferieur strict 


>= 


superieur ou egal 


<= 


inferieur ou egal 


.and. 


et 


.or. 


ou 


.not. 


non 



- les variables de type logical peuvent aussi etre uti- 
lisees 

- le else if et le else sont optionnels 
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f si x<0 ou x>2 

autre exemple : calcul de y = < y/x si xg[o,i] 

^1 si x>l 



rac:if (x<0.0. or .x>2.0) then 

y = 0.0 
elseif (x< = 1.0) then 

y=sqrt (x) 
else 

y=i.o 

end if rac 
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2.7.2 La selection select case 

- pour choisir une instruction a executer en fonction 
de la valeur d'une expression scalaire de type entier, 
caractere ou logique 

- exemple : 

cha racter ( len =30) pays 
langue : se I e c t case(pays) 

case( 'france ' , ' quebec ' , 'suisse ' , 'belgique ') 

print*, ' bonjour ' 
case ( ' royaume— uni ' , ' us a ' ) 

print*, 'hello ' 
case default 

p ri n t *,' langue pas disponible ' 
end select langue 



- l'expression testee doit verifier au plus un case 

- le case default n'est pris en compte que si l'expres- 
sion ne verifie aucun case 
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2.7.3 Iterations : structure do 



forme enumerative (nombre d'iterations fixe) 



n 



exemple : calcul de y = 



x 



1=1 



integer 
real x , y 



1 ,n 



y=o.o 

do i = l,n 

y=y+x * * i / i 
end do 



syntaxe 



nom 



do var iable=debut , f in [ , pas 
instructions 



end do 



nom 



autre exemple : calcul de 5 = 1 + 3 + 5 + 7 + 9 



integer 



1 , s 



s=0 

do i = l,9,2 

s=s+i 
end do 
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forme infinie : on ne sort qu'avec l'instruction exit 



nom 



end do 



do 

instructions 

if (condition) then 

exit [nom] 
end i f 

instructions 
nom 



exemple 



do 

read * , nombre 

if (nombre==0) then 

exit 
end i f 

somme = somme + nombre 
end do 
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- forme conditionnelle : arret si la condition specifiee 
en debut de boucle est satisfaite 



c 


nom : 
snd d( 


c 


io w 

Lnsti 
nom 


hile (condition) 
ructions 


exemple : 



reponse= ' oui ' 

lecture : do wh i le ( reponse= ' oui ' ) 



pri n t * , ' continuer ?' 
read * , reponse 
end do lecture 
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2.8 Les entrees-sorties ecran/clavier 

lecture au clavier de la valeur d'une variable, quel 
que soit son type 
read*, variable 

ecriture a l'ecran d'un message et/ou de la valeur 
d'une variable 
print*, 'bonjour' 
print *,x 

print*, ' il est ' ,h, 'heures' 
- exemple : 



integer n 

character ( len =15) pays 

print*, 'quel age avez— vous ?' 
read * ,n 

print*, 'vous avez ',n,'ans' 
print*, 'de quel pays etes vous ?' 
read * , pays 

print*, 'vous venez de ' ,pays 
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2.9 Compilation 

- pour faire fonctionner un programme ecrit en f90, il 
faut le traduire en un langage comprehensible par 
l'ordinateur : c'est la compilation. 

cette compilation est effectuee par un compilateur. 
Dans ce cours, nous utiliserons le compilateur libre 
gfortran. 

- la compilation cree un fichier dit executable. 

- exemple : le programme s'appelle calcul, il est ecrit 
dans le fichier toto.f90. On souhaite compiler ce pro- 
gramme pour creer un executable nomme ici com- 
pute. On tape done la commande 

gfortran -o compute toto.f90 

- plus generalement, la commande de compilation est : 

gfortran -o executable fichier-fortran 
ou executable est le nom que Ton souhaite donner au 
fichier executable et fichierJortran est le nom du fi- 
chier qui contient le programme f90 (il est conseille 
de donner au fichier le meme nom que le programme, 
en rajoutant l'extension .f90). 

- differentes options existent, cf. section 3. 
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- messages d'erreurs : 

— ► le compilateur detecte certaines erreurs (syntaxe, 
conformite de types, de dimensions) 

—> il affiche un message donnant le numero de la ligne 
contenant l'erreur, avec une breve explication en 
anglais 

— > exemple . (programme contenu dans le fichier erreur_exemple.f90) 

erreur_exemple . f 90 : 7 . 3 : 

y = x*log(x) 
1 

Error: Symbol 'y' at (1) has no IMPLICIT type 

— ► une erreur en provoque souvent d'autres (exemple : 
une variable non declaree va generer des erreurs a 
chaque ligne ou elle est utilisee) 

— > la liste des erreurs est done souvent tres longue, et 
il faut alors commencer par corriger les premieres 
de la liste, puis recompiler. 
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> les messages d'erreurs ne sont pas toujours tres 
clairs ... 

exemple : dans un programme contenu dans le fi- 
chier erreur_declaration.f90, on fait une erreur de 
frappe dans la declaration suivante : 

ineger : : n 
au lieu de : 

integer : : n 
ce qui provoque le message suivant : 

erreur_declaration . f 90 : 5 . 2 : 

ineger : : n 
1 

Error: Unclassif iable statement at (1) 
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2.10 Execution 

- si le fichier contenant l'executable s'appelle 



prog 



alors taper simplement ./prog 



erreurs a l'execution 

—> toutes les erreurs ne sont pas detectees a la com- 
pilation (certaines options peuvent cependant y 
aider) 

— > divisions par 0, erreurs d' affect at ion de tableau 

— ► erreurs d'algorithme (+ au lieu de — ) 

— >• il faudra apprendre a trouver ces erreurs en fai- 
sant des tests, et en utilisant un debugger (cf. 
sec. 3.5.5). 
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2.11 Ecrire un programme proprement 

- ecrire des commentaires (but et fonctionnement du 
programme, signification des variables, decrire chaque 
bloc d'instruction) 

- indenter correctement (en particulier les structures 
de controle) 

- nommer les boucles imbriquees 

- declarer les variables dans l'ordre croissant des struc- 
tures (logiques, entiers, caracteres, reels, puis tableaux) 
mieux vaut utiliser des noms de variables longs mais 
clairs que des noms courts et obscurs 

- il est fortement conseille d'utiliser l'editeur de texte 
Emacs pour programmer. Outre sa puissance d 'edition 
en general (au moyen de raccourcis clavier qui per- 
mettent de taper du texte tres rapidement), il possede 
des fonctions adaptees a la programmation : inden- 
tation automatique, completion des noms de pro- 
grammes et de boucles, coloration des mots clefs, 
etc. Voir le guide Emacs distribue en cours [5]. 
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2.12 Exemple complet 


Till 1 j • i — x n l 'xl l l 

calcul de la solution de x — e = par la metnode de 


Newton. 




algorithme : 




xq et tol donnes 




calculer x n+1 = x n 

1 + e Xn 




tant que > tol 




program 


zero 






i m p 1 i c 


it none 






/ variables 






real 


: x , y , dif f , tol 






/ lecture des donnees 






print* 


,'entrer le x initial' 






read * , 


X 






print* 


, ' t o 1 erance ' 






read * , 


tol 






/ methode de newt on 






newt on 


: do 






/ -- 


i t erat ion 






y=x 


-(x-exp(-x))/(1.0+exp(-x)) 






/ -- 


difference entre deux iteres 
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if (x/ = 0.0) then 

diff=abs ((x-y)/x) 

else 

dif f=abs (x— y) 

end i f 

/-- mise a jour 
x=y 

/-- test de convergence 
if (diff<tol) then 

exit 
end i f 

end do newton 

print *,' solution approchee= ' , x 
end program zero 
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3 Un peu plus loin ... 

3.1 Variables "const antes" 

Pour declarer une variable "constante" (c.-a-d. qui a 
une valeur non modifiable) : 

- attribut parameter 

- la variable est non modifiable, sinon erreur detectee 
a la compilation 

- utile pour declarer et proteger des constantes phy- 
siques 

- exemple : 

real, parameter :: pi=3. 14159265, N=6.02252e23 

dans la declaration, les operations usuelles sont uti- 
lisables, mais pas l'appel de fonction : 
real, parameter :: quart=l. 0/4.0 / — vaiide 
real, parameter :: pi=acos(— 1.0) / — non vaiide 
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3 . 2 Sous-programmes 

3.2.1 Objectif 

- bloc d'instructions utilise plusieurs fois =4> l'ecrire 
dans un sous-programme une fois pour toutes, et 
l'appeler quand necessaire (aussi appele procedure) 
deux types de sous-programme : les subroutines et 
les f onctions 
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— exemple 1 : subroutine 




program exproc 






implicit none 






■ 

real . . moyenne , maximum 
real , dimension ( 100) tab 






call random number(tab) 






call sp ( tab , moyenne , maximum ) 






print *, moyenne , maximum 






contains 






subroutine sp( t, moy, max ) 






implicit none 






real , dimension ( 1 00) , intent(in) t 






i ■ ■ ■ / ■ \ 

real , intent(out) moy, max 






i n teger i 






max — t(IJ; moy = t(lj 
do i = 2,100 






if (t(i) > max) then 






max — t ( i ) 






end i f 






moy = moy + t ( i ) 
end do 






mov = mov/100 






end subroutine sp 






end program exproc 
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— exemple 2 : fonction 




program exfct 






implicit none 






real : : maximum 

real , dimension ( 100) tab 






call random number(tab) 






maximum = maxi(tab) 
print * , maximum 






contains 






function maxi (t ) 






implicit none 






real , dimension ( 1 00) , intent(in) t 
integer i 

real :: maxi /- la fonction est declaree 






maxi = t ( 1 ) 
do i = 2,100 

if ( t(i) > maxi ) then 

maxi = t ( i ) 
en d if 
end do 
end function maxi 






end program exfct 
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3.2.2 Fonctionnement 

- place dans le programme : 

avant la derniere ligne end program nom_programme 
et apres le mot-clef contains 

program nom_programme 
declarations 
instructions 
contains 

subroutine nom_subrout ine 

end subroutine nom_subrout ine 

function nom_f onct ion 

end function nom_f onct ion 
end program nom_programme 

on parle alors de sous-programme interne (par oppo- 
sition aux sous-programmes externes, cf. section 4.3) 

- portee des variables : toutes les variables du pro- 
gramme principal sont connues de ses sous-programmes 
internes. Pas besoin de redeclaration, sauf pour les 
arguments (voir plus loin). 

les variables utilisees seulement par un sous-programme 
(variables locales) : declarees a l'interieur de celui-ci, 
comme dans un programme principal. Connues de 
lui seul. 
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- appel d'un sous-programme : 

— > subroutine : appelee par l'instruction call 

call nom_subrout ine ( argument s ) 

— > fonction : appelee comme une fonction mathe- 
matique : 

y=nom_f onct ion ( arguments ) 



3.2.3 Arguments d'une subroutine 

subroutine : sous-programme ayant ou non des argu- 
ments d'entree et de sortie : 

- variables a declarer a l'interieur de la subroutine 

- il est conseille de preciser dans ces declarations les 
attributs suivants 

intent (in) pour les arguments d'entree 
intent (out) pour les arguments de sortie 
intent (inout) pour les arguments mixtes. Ainsi le 
compilateur detecte une mauvaise utilisation des ar- 
guments (entree modifiee, sortie non donnee, etc.) 
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3.2.4 Arguments d'une fonction - declaration 
du resultat 

fonction : comme une subroutine, mais 

- arguments d 'entree obligatoires 

renvoie necessairement un resultat, stocke dans une 
variable ayant le meme nom que la fonction 
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- le type du resultat peut etre donne soit avant le nom 
de la fonction (sauf si resultat tableau) : 




real function maxi(t) 






end function maxi 




soit dans une declaration : 






function maxi (t ) 






real maxi 






end function maxi 




- le resultat peut avoir un nom different du nom de la 
fonction, en utilisant la clause result 




function maxi(t) result(y) 






real y 






end function maxi 




dans ce cas, pas de type avant le mot clef function. 
- la variable resultat ne doit pas avoir d'attribut intent 
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3.2.5 Les tableaux en arguments 

probleme : dans l'exemple de la page 34, la subrou- 
tine sp n'accepte que des tableaux de taille 100. 
Comment faire pour qu'elle fonctionne avec des ta- 
bleaux de tailles quelconques ? 

le plus simple : donner en entree le tableau et ses 
dimensions ("assumed size array" en anglais) : 



program exproc 




real , dimension ( 100) 


: : "t ab 


rnntr)in<! 

\_- \J II t U 1 II J 




subroutine sp( t, n , 


moy , max ) 


integer , intent ( in 


) :: n 


real , dimension ( 1 : 


n) , intent ( in ) t 


real , intent ( out ) 


: : moy , max 


integer i 




max — t ( 1 ) ; moy — 


t(i) 


do i=2,n 




end do 




moy = moy / 100 




end subroutine sp 




end program exproc 
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autre exemple 



program toto 

real, dimension ( 1 : 3 , 1 : 2 ) A 

call sousprog ( A , 3 , 2 ) 

contains 

subroutine sousprog (A , m , n) 

integer, intent(in) m , n 
real , dimension ( 1 : m , 1 : n) A 

end subroutine sousprog 
end program toto 
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- on peut aussi ne passer que le tableau en argument 
et recuperer ensuite ses dimensions (tableau a profil 
implicite, " assumed shape array" en anglais) . Atten- 
tion alors a la declaration : 




program exproc 






real , dimension ( 100) tab 






contains 

subroutine sp( t, moy, max ) 






real, dimension (:), intent(in) t 
real , intent out j moy, max 
integer i , n 






n=size (t ) 

max — t ( 1 ) ; moy — t ( 1 ) 
do i=2,n 






end do 

moy = moy / 100 






end subroutine sp 






end program exproc 
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- autre exemple : 



program toto 




real , dimension ( 1 : 3 , 1 : 2 ) 


: A ,B 


real , dimension ( 1 : 3 ) u, 


V 


/** *~\ 1 /~\ n vn "v** /-\ / A T2 i 

cdii sousprog 1 ii , D ) 




V — I C t ^ U J 




contains 




subroutine sousprog ( A , B ) 




real , dimension ( : , : ) 


, intent ( in ) A 


real , dimension ( : , : ) 


, intent ( out ) B 


integer m,n 




m= s i z e y A , 1 ) ; n= size 


(A ,2) 


B=A * * 2 




end subroutine sousprog 




function fct(u) result(v 


) 


real , dimension ( : ) , 


intent ( in ) u 


real , dimension ( size 


(u)) :: v 


end function fct 




end program toto 





ATTENTION : le profil implicite s'utilise seulement 
pour les tableaux arguments. Pour des tableaux lo- 
caux ou un resultat de fonction dont les tailles de- 
pendent de celles des arguments, remarquer l'utili- 
sation de l'instruction size. 
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3.2.6 Les chaines en arguments 

meme principe que pour les tableaux : 

- on peut passer en argument la chaine et sa longueur 

- on peut aussi utiliser un profil implicite, pour les 
arguments uniquement 

- exemple : 

program chaine 

character ( len =20) :: ch 

print * , ch , f 2 ( ch ,20) 
print*, f 1 ( ch) 

contains 

function f 1 ( ch) 
implicit none 

character ( len =*) , intent(in) ch 
character ( len = len ( ch )) ) :: fl 

end function eff_blanc 

function f 2 ( ch , n) 
implicit none 
integer , intent(in) n 
character ( len=n) , intent(in) ch 
integer f 2 

end function nb blancs 



end program chaine 
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3.3 Tableaux 

3.3.1 Tableaux statiques 

deja vu : les dimensions sont fixees a la declaration : 
real, dimension(l:10,l:4) :: A 

3.3.2 Allocation dynamique 

- il est frequent que les tailles des tableaux utilises 
dans un programme dependent de donnees et ne 
soient pas connues a la compilation 

- moyen simple et efficace : attribut allocatable 
(tableau allouable) : 

real , dimension ( : ,:) , allocatable :: A 

- le nombre de dimensions est donne (ici 2), mais pas 
les dimensions elles-memes 

- pour manipuler ces tableaux, on commence par les 
allouer en memoire (c.-a-d. la place necessaire est 
reservee), en donnant leurs dimensions. Exemple : 

integer n , m 
read * ,n , m 

a I locate ( A ( 1 : n , 1 : m) ) 
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- pour changer la taille du tableau : le desallouer, puis 
le reallouer : 

deallocate (A) 

allocate (A(l:3 ,0:4)) 

- evidemment, on ne peut utiliser un tableau non en- 
core alloue. 

- ne pas oublier de liberer la place memoire avec deallocate 
quand le tableau n'est plus utilise (en particulier 
pour les tableaux locaux dans des sous-programmes, 
sinon, risque de depassement de memoire). 

3.3.3 Allocation dynamique dans un 

sous- programme : tableaux automatiques 

- autre moyen simple de creer un tableau dont la taille 
n'est pas connue a la compilation 

- utilisable seulement au sein d'un sous-programme 

- principe : les dimensions du tableau sont des va- 
riables passees en arguments du sous-programme ("au- 
tomatic array" en anglais). Exemple : 

subroutine autom(m,n) 

integer m , n 

real , dimension (m , n) A 



46 



- avantage : plus rapide en execution que la creation 
d'un tableau allouable. A utiliser pour les sous-program- 
mes appeles tres frequemment. Les instructions for- 
tran sont aussi plus simples a ecrire. 

- inconvenient : aucun moyen de verifier s'il y a assez 
de place pour creer le tableau, contrairement aux 
tableaux allouables. Le message d'erreur correspon- 
dant est generalement tres obscur. 

- pour les tableaux de grandes tailles utilises peu frequem- 
ment, il est done plutot conseille d'utiliser les ta- 
bleaux allouables. 

3.3.4 Terminologie 

- rang : nombre de dimensions (le rang de A est 2) 
taille d'une dimension : nombre d'elements dans cette 
dimension (les tailles des dimensions de A sont 3 et 

5) 

forme : suite ordonnee des tailles des dimensions (la 
forme de A est 3,5) 

- une taille peut etre nulle 



47 



3.3.5 Operation conformes entre tableaux 

but : faire des operations sur des tableaux sans ecrire 
les operations element par element 

- deux tableaux sont conformes s'ils ont meme forme : 
A (1:10,1:4) et B (3:12,0:3) sont conformes 

A (1:10,1:4) et B (3:12,1:3) ne sont pas conformes 
A (1:1,1:1) et B(l:l) ne sont pas conformes 

- une expression est conforme si elle fait intervenir des 
tableaux conformes 

- un scalaire est considere comme conforme a tout ta- 
bleau 

la plupart des operations scalaires sont applicables 
aux tableaux : operation appliquee element par element, 
resultat de meme forme que les tableaux operandes 
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- exemples : 



program conforme 




implicit none 




integer 1 , j 




real , dimension (3 ,3) 


: A , B , C 


call random.number (A) 




II i i / rt \ 

call random.number (Bj 




C=A+B + 2.0 /-- op. 


equival ent e a 


print * ,C 




do i = l,3 




do j =1,3 




C(i, j)=A(i, j)+B( 


i, j) + 2.0 


end do 




end do 




print * ,C 




C=A*B /-- pas 


prod, matriczel 


p r i n t * , C 




C=sqrt (B)*C 




p r i n t * , C 




C=B+cos (B*3.1416) 




p r i n t * , C 




A=l-sin( loff(B+A)) 




p r i n t * , A 




end program conforme 
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3.3.6 Creer des tableaux "a la main" 





avec une liste de variables d'un meme type 

varl var2 varn on dpi it creer nn vecteur ft a- 

bleau de rang 1) avec le construct eur (/.../) : 


V 


= (/varl , var2 , . . . , varn/) 




ia nsie peui-ene ineiiie ene cieee avec une uu(iLit< 
implicite 


V 
V 


= (/ (i**2,i = l,4)/) .'--equivaut a 

= (/l,4,9,16/) 




on peut imbriquer plusieurs boucles implicites 


V 
V 


= (/((i+j , j =1 ,2) , i = l ,3)/) ! -equivaut a 

= (/2,3,3,4,4,5/) 
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- comme dans une boucle do on peut donner un pas : 

V = (/ (i,i = l,ll,2)/) /-- equiv aut a 

v = (/l,3,5,7,9,ll/) 

- pour transformer un vecteur en tableau a plusieurs 
dimensions : 

—> fonction reshape(vecteur,f orme_du_tableau) 

real, dimension (2 ,3) A 

A=reshape(v,(/2 ,3/)) 

donne la matrice A = (s7ii) ( r enrphe colonne par 
colonne) 

3.3.7 Les sections de tableaux 



on accede a une section de tableau en modifiant conve- 
nablement les indices de debut et de fin des dimen- 
sions : 



real . 


dimension ( 1 : 


10,1:20) :: A 


A(2:6 , 


1:20) /-- 


l i gn es de 2 d 6 


A (2:6 , 


:) 


idem 


A(i,:) 




i erne I i gn e 


A(: , j) 




j erne colonne 


A(2:6 , 


1:4) /- 


element communs aux 






I i gn es 2 d 6 et 






c o I o nn e s 1 d ^ 
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3.3.8 Stockage des tableaux dans la memoire 
et ordre des boucles 

- il est tres important lorsqu'on manipule de gros ta- 
bleaux de savoir comment ils sont stockes en memoire 

- exemple : y-a-t'il une difference importante entre les 
deux blocs suivants ? 

do i = l,m 

do j=l,n 

c(i, j)=a(i, j)+b(i, j) 

enddo 
enddo 



do j=l,n 

do i = l,m 

c(i, j)=a(i, j)+b(i, j) 

enddo 
enddo 



=4> oui : le premier bloc est beaucoup plus couteux 
en temps ! 
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quel que soit le nombre de dimensions d'un tableau, 
il est stocke en memoire sous forme uni-dimensionnelle. 
Les elements sont stockes ainsi : 



a(l,l) 


a(2,l) 


a(3,l) 




a(m, 1) 


a(l,2) 


a(2,2) 


a(3,2) 










a(l,n) 


a(2,n) 


a(3,n) 




a(m,n) 





exemple : les elements du tableau suivant 
integer, dimension(2,4) :: a 

sont ranges ainsi : 

a(l,l) a(2,l) a(l,2) a(2,2) a(l,3) a(2,3) a(l,4) a(2,4) 



a(l, 1) 


a(l,2) 


a (1, 3) 


a(l,4) 


a (2,1) 


a (2, 2 


a (2,3) 


a (2,4) 



a (1 


,D 


a (2 


, 1) 





a(l 


,4) 


a (2 


,4) 

f 



conclusion : quand on parcourt les elements d'un ta- 
bleau, il faut que la boucle sur la premiere dimension 
soit la boucle la plus interne (ceci est lie a la notion 
de cache dans la memoire) 

c'est l'inverse de ce qui se fait en C, et de ce qui se 
fait en general quand on ecrit les algorithmes 
voir en TP pour une comparaison en temps calcul 
des deux strategies. 
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3.4 Entrees-sorties (E/S) 

On a vu dans la premiere partie comment lire au clavier 
et ecrire a l'ecran. Ici, on va voir comment lire et ecrire 
dans un fichier. 



3.4.1 E/S en format libre, dans un fichier texte 
a acces sequentiel 

C'est le plus simple, equivalent de ce qu'on a vu avec 
l'ecran/clavier 

- ouverture du fichier toto : 
open (unit = 10,file = 'toto ' ) 



■> cette instruction ouvre le fichier a la premiere 
ligne (s'il n'existe pas, il est cree). 

■» l'expression unit=10 indique que le fichier toto 
est connect e au descripteur numero 10. 

■» tous les nombres peuvent etre utilises comme 
descripteurs sauf 5 (clavier) et 6 (ecran). 

si le mot clef file n'est pas donne, le fichier 
est cree, ou n est le descripteur donne 



fort.n 



apres le mot clef unit. 
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ecriture de donnees dans toto (si associe au descrip- 
teur 10) : 

write (10,*) donneel , donnee2 , donnee3 



> une instruction write ecrit toutes les donnees sur 
la meme ligne, le write suivant passera a la ligne 
suivante. 

> les donnees peuvent etre de n'importe quel type 
(constantes ou variables) 

> le print*, pour ecrire a l'ecran est equivalent a 
write (6,*) et a write (*,*) 
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- lecture de donnees dans toto (si associe au descrip- 
teur 10) : 

read (10,*) varl , var2 , var3 

—> la premiere ligne de toto est lue, le read suivant 
passera a la ligne suivante. 

— >• les donnees sont supposees etre separees par des 
blancs. 

—> elles sont stockees respectivement dans varl, var2, 
varn. 

—> si plus de donnees que de variables : seules les 
premieres donnees sont stockees (cas inverse : seules 
les premieres variables sont affectees) 

—> le read*, pour lire au clavier est equivalent a 
read(5,*) et a read(*,*). 

- fermeture du fichier toto : 
close (10) 

—> attention, si on reouvre le fichier avec open, on se 
retrouve a la premiere ligne. 

- ecrasement de fichier : un write apres l'ouverture 
d'un fichier existant detruit le contenu du fichier. 



56 



- exemple 1 : 

fichier titi . txt : 


UOI1JUU.X 4: lie UI CO 

temperature 37.2 


programme : 


character ( len =1 4) c 1 c 2 c3 
integer n 
real : T 




open(unit=10,file = 'titi .txt ') 

I'D 3 n 1 1 D -4s \ r* i *n r* O / 1 o ot or* & 

ItrclUl-LUj^J C 1 , 11 , C Z ; L (z. L L Lb f C 

read(10,*) c3 , T /-- 
close (10) 


± kL ! H. u U ft C 

2eme ligne 


open(unit=20,file = 'toto.txt ') 




/-- ecriture 1 ere ligne 
write(20,*) cl , ' i 1 est n,c2 




/-- ecriture 2 erne ligne 
write(20,*) 'la ',c3,'est de ',T, 


' degres ' 


close (20) 





apres execution, fichier toto.txt : 



bonjour il est 4 heures 

la temperature est de 37.20000 degres 
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- exemple 2 : dans le fi drier info.txt, liste de 1000 per- 
sonnes avec un nom, age et taille sur chaque ligne. 
Programme de lecture : 

integer i 

cha racter ( len =30) , dimension ( 1 : 1 000 ) nom 
integer, dimension ( 1 : 1 000 ) age 
real , dimension ( 1 : 1 000 ) taille 

open(unit=30,file = 'info .txt ') 
do i = l,1000 

read(30,*) nom ( i ) , age(i), taille(i) 
end do 
close (30) 



58 



3.4.2 E/S dans un fichier binaire a acces 
sequentiel 

les fichier s de textes sont des fichier s de caracteres, ce 
qui peut prendre beaucoup de place s'il y a beaucoup 
de donnees. 

- dans ce cas, il vaut mieux stocker les donnees (en 
particulier numeriques) sous forme binaire. Le fichier 
n'est pas lisible avec un editeur de texte, mais il est 
beaucoup plus petit. 

- ouverture : avec open comme pour un fichier texte 
en precisant le parametre form = 'unformatted' : 

open(unit=15,file = 'toto ' , form=' unformatted ') 

- lecture/ ecriture : open et write sans le parametre * : 

read (15) varl , var2 
write (15) donnee 
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- exemple : tableau de reels aleatoires, stocke dans un 
fichier texte et dans un fichier binaire. Comparer les 
tailles des deux fichiers. 



real , dimension ( 1 : 1000) :: tableau 

call random. number ( tableau ) 

open(unit=l,file = 'texte.txt ') 
write(l,*) tableau 
close ( 1) 

open(unit=2,file = ' binaire. b ' , form=' unformatted ') 
write (2) tableau 
close (2) 
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3.4.3 E/S formattees 

- inconvenient de E/S en format libre : 

— > fichier pas propre : reels ecrits sous forme decimale 
si > 10 -2 , sous forme scientifique sinon (cf le fi- 
chier texte.txt de la section 3.4.2). 

—> lecture peu sure : risque de melange des variables 

- format : on indique precisement le nombre et le type 
des variables a lire et ecrire 

- parametre fmt de read et write : 

read(descripteur , fmt= '(chaine) ') var 1 , . . . 
write (descripteur , fmt= '(chaine) ') var 1 , . . . 

ou chaine indique le type des donnees dans l'ordre 

ou elles sont lues ou ecrites 

- formats de base (plus de details sur [3]) : 

— > fmt='(nlm)' : n entiers de m chiffres 

—> fmt='(nFm.d)' : n reels sous forme decimale, formes 
de m caracteres (dont le .), dont d apres la virgule 

—> fmt='(nEm.d)' : n reels sous forme scientifique, formes 
de m caracteres, dont d apres la virgule (dont le .) 

—> fmt='(nAm)' : n chaines de m caracteres 

Dans tous les cas, c'est le caractere espace qui separe 
les n champs. 

- exemple : 
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character ( len =10) pr 


enom , nom 


integer age 




prenom='JEAN ' 




nom= 'DURAND' 




age=20 




write (* ,fmt=' (2A10,1 13) 


' ) prenom , nom , age 


va ecrire : 


JEAN__DURAND_. 
ALAIN _._MARTIN__ 


...20 
...33 



- si on doit utiliser plusieurs fois un meme format, 
on peut stocker la chaine associee au parametre fmt 
dans une variable. Exemple : 



f ormate=' (2A10,1 13) ' 

write (*, fmt=formate ) prenom, nom, age 
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3.4.4 Controle des erreurs en E/S 

avec les instructions precedentes, le programme s' arret 
a la moindre erreur (lecture dans fichier inexistant, 
lecture depassant la taille du fichier, ecriture texte 
dans un fichier binaire, etc.) 

des parametres des commandes open, read et write 

permettent d'eviter l'arret et de determiner la source 
de l'erreur. 
- exemple : le parametre iostat : 



integer erreur 

read ( 15 , * , iostat=erreur ) /- fich. texte 
read (23 , iostat=erreur ) / - fich. binaire 



-> si lecture correcte =4> erreur vaut 

■» sinon, si erreur < alors fin de fichier ren- 



contree, sinon autre erreur. 
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— > exemple : 
erreur=0 

do while ( erreur==0) 

read ( 15 , * , iostat=erreur ) var 
print*, 'variable=' ,var 

end do 

if (erreur<0) then 

print*, 'fin du fichier ' 
else 

pri nt * , ' erreur de lecture' 
end i f 



> voir les instructions open et write pour la signi- 
fication de iostat. 
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3.5 Trouver des erreurs dans un pro- 
gramme 

3.5.1 Les eviter et les detecter 

1. avant de programmer un algorithme, s'assurer qu'il est 
correct 

2. programmer proprement (voir sec. 2.11), utiliser au 
maximum les signaux d'erreurs renvoyes par les com- 
mandes (parametre iostat, etc.), afficher des messages 
d'erreur en cas d'operation interdite (y/x < 0, etc.) 

3. lors de la phase de developpement, utiliser des options 
de compilation pour detecter les erreurs a l'execution 

qui detecte en particulier les de- 

qui 



comme 



-fcheck=all 



passements de tableau) ou -ffpe-trap=invalid,zero 



detecte les operations invalides 

4. effectuer des tests judicieux pour tester le programme 

5. en desespoir de cause, raj outer des affichages de va- 
riables dans le programme, et utiliser l'instruction stop 
pour que le programme s'arrete a l'endroit voulu. 

La tache 5 peut etre simplifiee en utilisant un "debugger" 
(voir sec. 3.5.5). 
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3.5.2 Erreurs detectees a la compilation 

le compilateur renvoie un message plus ou moins ex- 
plicite (voir sec. 2.9) 
- erreur de syntaxe : exemple la compilation sur natty 
de 

program erreur_syntaxe 
implicit none 

integer : x 
real r 

r = l(T2 

end program erreur.syntaxe 

renvoie les messages suivants : 

erreur_syntaxe . f 90 : 4. 11 : 

integer : x 
1 

Error: Invalid character in name at (1) 
erreur_syntaxe . f 90 : 7 . 2 : 

r=10~2 
1 

Error: Unclassif iable statement at (1) 
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- oubli de declaration : voir l'exemple precedent 

- erreur de dimension : la compilation sur natty de 

program erreur_dimension 
implicit none 

real , dimension ( 1 : 4) a 

a(4,4) = 1.0 
end program erreur_dimension 



renvoie les messages suivants 

erreur_dimension . f 90 : 6 . 3 : 

a(4,4)=1.0 
1 

Error: Rank mismatch in array reference at (1) (2/1) 

- malheureusement, il existe beaucoup d'autres erreurs 
possibles ... 
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3.5.3 Erreurs detectees a l'execution 

II y a evidemment beaucoup d'erreurs possibles. En 
voici quelques unes avec le programme suivant 

program erreur_execut ion 
implicit none 
real x 

x=l/3 
print * ,x 

x=-1.0 

print * , sqrt (x) 
read * , x 

end program erreur_execut ion 

En supposant que l'utilisateur tape la lettre a au cla- 
vier, on aura alors les amchages suivants a l'ecran : 

0.0000000 

NaN 

a 

At line 11 of file erreur_execution.f 90 (unit = 5, file = 'stdin') 
Fortran runtime error: Bad real number in item 1 of list input 
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- le premier print renvoie 0, car la division est entiere 
(oubli de .0 sur le numerateur et/ou denominateur) . 
Mais il n'y a pas d'erreur d'execution a proprement 
parler, puisque le programme continue a tourner) 

- le deuxieme renvoie NaN (not a number), car la racine 
d'un nombre negatif n'est pas definie (la non plus, 
pas vraiment de message d'erreur : le programme 
continue a tourner) 

le dernier amchage est un message d'erreur : le pro- 
gramme s'arrete (on dit qu'il "plante") a cause de la 
lecture d'un caractere au lieu du reel attendu. 

3.5.4 Erreurs d'algorithme 

Ce sont des erreurs souvent dimciles a trouver. C'est 
pourquoi il vaut mieux bien ecrire son algorithme avant 
de le programmer 
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3.5.5 Les "debuggers" 

La tache 4 decrite a la section 3.5.1 peut etre simplifiee 
avec les logiciels dits "debuggers" : 

- sans eux, il faut recompiler et executer le programme 
apres chaque modification, ce qui peut etre long quand 
le programme est gros ; 

- un debugger permet, sans modifier le programme, 
d'afHcher la valeur des variables au cours de l'execution, 
et de stopper temporairement l'execution a l'endroit 
desire. 

l'utilisation d'un debugger necessite un certain inves- 
tissement, mais cet effort est largement recompense 
par la suite. 

Sur la machine natty, le debugger est gdb. II peut etre 
utilise avec l'interface graphique ddd (voir le manuel [2] 
avec un exemple simple et instructif dans la section 
"sample session"). 



Important : pour utiliser un debugger, il faut compiler 
le programme avec l'option I -g 
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4 Programmation modulaire 

4.1 Introduction 

4.1.1 Qu'est-ce-que c'est ? 

- programmer de facon modulaire consiste a decouper 
le plus possible un programme en morceaux algorith- 
miquement independants (appeles sous-programmes) 

- exemple : calculer t = tiSice(A~ 1 B). Ce probleme 
contient trois calculs qui peuvent etre programmes 
independamment : 

A i — ► trace(A) 
(A, B) i — y AB 
A^A' 1 

on essaiera done de rendre ces morceaux les plus 
independants possible au niveau informatique. 

4.1.2 Interet 

- permet de reutiliser plusieurs fois un morceau de pro- 
gramme facilement 

- le programme est plus lisible 

- le risque d'erreurs de programmation est moins grand : 
chaque bloc est ecrit, compile, et teste separement. 
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- exemple : creation d'un ensemble de fonctions pour 
faire des calculs matriciels (produit mat rice- vecteur, 
triangularisation, solution d'un systeme lineaire, de- 
terminant, trace, etc.) 

on a deja vu un moyen de creer des sous-programmes 
en utilisant les fonctions et subroutines. On va aussi 
voir les modules. 

4.2 Les procedures internes 

- deja vu a la section 3.2, facile a faire 

- probleme : pas emboitable. Une procedure interne 
ne peut pas contenir de sous-programme. Or une 
procedure interne peut necessiter l'utilisation d'une 
autre procedure interne (exemple : la procedure de 
resolution d'un systeme lineaire par methode de Gauss 
necessite une procedure interne de resolution d'un 
systeme triangulaire) 

- probleme de lisibilite : le programme principal peut 
etre tres gros car il contient touts les sous-program- 
mes 

- cette methode ne repond done pas parfaitement au 
probleme. 
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4.3 Les procedures externes 

4.3.1 Pourquoi "externes?" 

Ces procedures sont ecrites a P exterieur du programme 
principal : 

- soit apres 1'instruction end program qui termine le 
programme principal 

- soit dans un autre fichier 

4.3.2 Utilisation : probleme d'interface 

Elles peuvent s'utiliser exactement comme une procedure 
interne (meme syntaxe), mais il existe un probleme 
d'interface : 

- programme principal et procedures sont compiles se- 
parement 

- il n'y a pas de controle de coherence des arguments 
entre Vappel et la procedure appelee : les erreurs 
eventuelles ne sont done pas detectees a la compi- 
lation. Elles seront "visibles" a l'execution, mais se- 
ront peut etre difficiles a localiser. Exemple : types 
des variables differents entre l'appel et la procedure 
appelee. 
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- on dit que V interface entre programme principal et 
procedure appelee est implicite (contrairement aux 
procedures internes ou l'interface est dite explicite). 

- consequence : il est impossible de passer des tableaux 
ou des chaines a profils implicites dans une procedure 
externe avec interface implicite. 

=4> pour rendre une interface explicite, le plus simple 
est d'enfermer la procedure externe dans un module 
(il existe d'autres solutions moins simples avec les 
bloc interface, voir [3]). 
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4.4 Les modules 

- utilises pour que l'interface entre programme princi- 
pal et procedure externe soit explicite 

- on enferme la procedure externe dans un module de 
la facon suivante : 






module nom_du_module 






implicit none 






contains 






subroutine nom_subrout ine 






end subroutine nom_subrout ine 






end module nom_du_module 




- le module est alors place dans un fichier (different 
de celui du programme principal) qui a pour nom 
nom_du-module.f90 

- dans le programme utilisant la subroutine, il suffit 
de rajouter l'instruction use nom_du_module en tout 
debut de programme, avant Tinstruction implicit none. 
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- exemple : erreur d'argument dans l'appel d'une pro- 
ceaure externe 


1. interface implicite (sans module) 
fichier testl.f90 






program testl 






implicit none 






call somme ( ' oui ' , ' non ' ) 






end program testl 




fichier somme.f90 






subroutine somme(x,y) 






implicit none 






real , intent ( in ) x , y 
real s 






s=x+y 

print*, ' somme= ' , s 






end subroutine somme 




en entree : 2 chaines au lieu de 2 reels =4> erreur pas 
detectee a la compilation. A l'execution, on a un 
resultat aberrant. 
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z. interface explicite avec module 


flchier test 2. ISO 




program test2 






use mod_somme 






implicit none 






call somme ( ' oui ' , ' non ' ) 






end program test2 




fichier modsomme.f90 




module mod_ somme 






implicit none 






contains 






subroutine somme(x,y) 






implicit none 






real , intent(in) x,y 
real s 






s=x+y 

print*, ' somme= ' , s 






end subroutine somme 






end module mod_somme 
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l'erreur est alors detectee a la compilation : 



test2.f90:7. 13: 

call somme ( ' oui ' , 'non' ) 
1 

Error: Type mismatch in argument 'x' at (1); passed CHARACTER (1) to REAL (4) 



78 



— autres utilisations aes modules 


—> permettent de definir des const antes, des variables, 
aes types aerives que i on peut utinser aans amerents 


sous-programmes sans les passer en argument. 
Exemple : 

Ti c ni Trt rtn f~t p> in n~tn tQD 
llLylllt;! 1 1 oULL—ULZs—l I LLLLI L.J c/U 




module mod_cte_math 






implicit none 






real, parameter :: pi=3. 141593 






end module mod_cte_math 




fichier aire. f 90 




program aire 






use mod_ct e_math 






implicit none 






real r , s 






r = 1.0 

s=p i * r * * 2 






end program aire 







79 



— ► permet aussi de passer des procedures en argu- 
ments, pour faire des fonctions de fonctions (voir 
section 5) 

- conseil : utiliser les modules autant que possible. 
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4.5 Compilation separee 

comment compiler un programme quand on a des sous- 
programmes externes situes dans des fichiers differents ? 

- compilation globale 

— > utiliser la commande vue sec 2.9 

gfortran -o nom_executable UsteJichiersJortmns 

—> mais les modules doivent etre compiles avant 
les sous-programmes utilisateurs. lis doivent done 
etre places avant dans la liste 

—> inconvenient : une modification d'un fichier en- 
traine la re-compilation de tous les fichiers 
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- compilation separee : 2 phases successives 

1. compilation (option -c) : pour tout fichier fortran 
(ici fi.ch.f90) faire 

gfortran -c fich.f90 

resultat : 

— ► creation du fichier objet fich.o (de meme nom 
que le fichier fortran mais d 'extension .o) 

—> si fich.f90 contient le module f ich, creation du 
fichier fich.mod (de meme nom que le module 
avec l'extension .mod). Contient l'interface du 
module, necessaire pour compiler les programmes 
utilisateurs de f ich 

— > attention : compiler les modules avant les sous- 
programmes utilisateurs 

2. edition de liens (option -o) : les fichiers objets (.o) 
sont lies ensemble pour faire l'executable 

gfortran -o nom_executable HsteJichiers-objets 
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- interet : 

— > on ne recompile que les fichiers modifies et les fi- 
chiers dependants, et Ton refait l'edition de liens 

—> outil complement aire efficace pour la compilation 
separee de gros programmes (plusieurs dizaines de 
fichiers) : la commande make (voir la fin de cette 
section) 
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Exemple : programme fortran90 compose d'un pro- 
gramme principal prog qui utilise deux modules modi 
et mod2. Le module modi utilise le module mod2. 



prog. 



(3) 



compilation 



pro; 



g.f90 



prog 



executable 



(4) 



edition de liens 



modl.o 

modl.mod 



(2) 



modl.f90 



fichiers 

mod2.o oblets 

mod2.mod 



(1) 



modlf90 



fichiers 
fortran 



=4> la compilation separee est : 

(1) compilation : 
f90 -c mod2.f90 

(2) compilation : 
{90 -c modLfOO 

(3) compilation : 
{90 -c prog. {90 

(4) edition de liens : 

{90 -o prog prog.o modl.o mod2.o 
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=> en cas de modification de prog, alors seules la 
compilation de prog.f90 puis l'edition de liens sont 
necessaires (phases (3) et (4)) 
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make 



: le fichier suivant 



utilisation de la commande 
(nomme "Makefile") contient les dependances entre 
fichiers et les actions a executer pour mettre a jour 
l'executable : 



# edition de liens 

prog : prog.o modl.o mod2.o 

<TAB> gfortran -o prog prog.o modl.o mod2.o 


# compilation 




prog.o : prog. f 90 


modi . o mod2 . o 


<TAB> gfortran -c 


prog. f 90 


modl.o : modi. f 90 


mod2 . o 


<TAB> gfortran -c 


modi. f 90 


mod2.o : mod2.f90 




<TAB> gfortran -c 


mod2.f90 



La commande unix 



make 



utilise ce fichier pour ne 



compiler que ce qui a ete modifie : il suffit de taper 
cette commande dans le dossier contenant le pro- 
gramme. 

Attention : la tabulation (visualisee ci-dessus par 
<TAB>) est indispensable. 
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- on peut parametrer ce "Makefile" pour utiliser des 
options de compilation differentes ou pour "nettoyer" 
le dossier : 




#parametres 
DEBUG= 

COMPILE=gfortran $ (DEBUG) 






# edition de liens 

prog : prog.o modl.o mod2.o 

<TAB> $ (COMPILE) -o prog prog.o modl.o mod2.o 
<TAB> Oecho comDilatioii terminee 






compilation 
prog.o : prog. f 90 modl.o mod2.o 
<TAB> $ (COMPILE) -c prog. f 90 
modl.o : modi. f 90 mod2.o 
<TAB> $ (COMPILE) -c modi. f 90 
mod2.o : mod2.f90 
<TAB> $ (COMPILE) -c mod2.f90 






# destruction des fichiers objets et modules 






detruire : 






<TAB> rm -f *.o *.mod 
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differents appels de make 



make 



: compilation normale 



make OPT=-g : compilation en mode debug 

: compilation optimise de niveau 2 



make OPT=-02 



make detruire 



: efface les .o et les .mod 



Voir [4, 1] pour des details. 



Remarque : pour les pet its programmes, le make- 
file peut etre genere a la main, mais pour les pro- 
grammes plus longs, il peut etre utile de le generer 
automatiquement, par exemple avec le script perl 
mkmf , voir : 



http: //www. gf dl .noaa.gov/ ~vb/mkmf .html 
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5 Utilisation avancee 

- les concepts presenters ici sont plus complexes que les 
precedents 

- a part le concept de precision des reels, les autres 
ne sont pas forcement necessaires pour realiser des 
programmes de calcul scientifique, sauf peut-etre les 
pointeurs qui peuvent etre indispensables dans cer- 
tains cas 

cependant, ces concepts peuvent permettre de sim- 
plifier l'ecriture des programmes 
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5.1 Types derives 

- structures de donnees renfermant plusieurs types 
differents 

- exemple : type personne contenant les nom, prenom, 
et age d'une personne : 

type personne 

character ( len =30) nom, prenom 

integer age 
end type personne 

on peut alors declarer des variables de ce type, comme 
pour n'importe quel type 

type ( personne ) :: etudiantl , etudiant2 
type (personne ) , dimension ( 1 : 3 5 ) td 
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- les different es variables d'un type derive sont ap- 
pelees champs. On y accede avec le caractere %, et 
on peut utiliser le constructeur de meme nom que le 
type pour les initialiser : 

etudiant l=per sonne ( 'Cesar' , 'Jules ' ,57) 
etudiant 2=etudiant 1 
print*, etudiant 2 

do i = l,2 

print*, ' nom , prenom , age ' 

read * , td ( i) %nom , td ( i) %prenom , td ( i)%age 
end do 

p r i n t * , td ( 1 ) 

l'affectation etudiant2=etudiantl recopie tous les 
champs. 
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5.2 Precision des reels 

- par defaut, le type real est dit "simple precision", 
c.-a-d. code sur 4 octets (=32 bits), avec 6 chiffres 
significatifs : ceci est souvent insuffisant pour faire 
des calcul scientifique precis ; 

il faut alors utiliser la representation "double precision" , 
c.-a-d. codee sur 8 octets (=64 bits), avec 15 chiffres 
significatifs => 3 possibilites : 

—> type double precision : 
double precision x 

—> type real avec parametre kind : 

integer, parameter pr=kind ( 1 . dO ) 

real (pr ) : : x 

(type reel de meme precision que l.dO, c.-a-d. 
double precision) 

— ► ou 

integer, parameter pr=s e I e c t e d _r e a I _k i n d ( 1 5 , 3 ) 

rea I (pr ) : : x 

(type reel avec 15 chiffres significatifs et exposant 
a 3 chiffres, c.-a-d. double precision) 

- la 3 e possibilite est la plus "portable" (donne la meme 
precision quelle que soit la machine) 
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- les valeurs numeriques s'ecrivent differemment selon 
la precision choisie 



valeur 


real 


double precision 


real (pr) 


1 

4,32 

4,32 x 10" 4 


1. 

4.32 

4.32e-4 


l.dO 

4.32d0 

4.32d-4 


l._pr 

4.32_pr 

4.32e-4_pr 



- attention a la coherence des types dans les affecta- 
tions ou operations : 



integer 


n 








double precision 




: u 


n = 1 










u = n 


/ — 


u 


= 1 


. 00000000000000 


u = 1.2d0 


/ — 


u 


= 1 


.20000000000000 


n = u 


/ — 


n 


= 1 




u = 1.2 


/ — 


u 


= 1 


. 20000004768372 



- par prudence : ne pas melanger les entiers et reels 

- complexes : comme pour les reels : 



complex z 




/ — 


si mp I e 


V ■ 


double complex 


z 


/ — 


double 


V ■ 


pr=kind (l.dO) 




/ — 


o u 




pr=selected _rea 1 _ 


kind (15 


,3) 






complex (pr ) z 




/ 


double 


V • 
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5.3 Fonctions de fonctions 

- il s'agit simplement de creer des procedures dont les 
arguments peuvent etre d'autres procedures 

- exemple simple : ecrire une fonction integrale qui 
pour toute fonction / et tout couple de reels a, b 
calcule une approximation de J a f(x) dx 

- pour cela : le fonction integrale doit etre mise dans 
un module, et le type de la fonction / et de ses argu- 
ments doivent etre declares dans un bloc interface 
a l'interieur de la fonction integrale 

- attention : dans le programme utilisant la fonction 
integrale, la fonction / passee en argument est : 

—> soit une fonction intrinseque (comme cos, log, 
exp, sqrt ...), il faut alors la declarer avec l'at- 
tribut intrinsic ; 

—> soit definie dans un module toto, il suffit alors 
de declarer l'utilisation de ce module (use toto) ; 

— > soit definie dans un sous-programme externe non 
place dans un module, il faut alors la declarer 
avec l'attribut external. 
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— exemple : 




module mod_ int egrale 
implicit none 






contains 






real function int egrale ( f , a , b ) 
implicit none 






/ argument s 

real , intent(in) a,b 






/ argument f : bloc interface 

interface 

rea 1 function f (x) 

real , intent ( in ) x 

end function f 
end interface 






/ var locales 

integer n , i 
real h 






n=100 ; h=(b-a)/n 
integrale =0. 
do i = 0,n—l 

int egrale=int egrale+h*f ( i *h) 
end do 






end function integrale 






end module mod_ int egrale 
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- utilisation de cette fonction dans un programme : 

program pr og_ int egrale 

use mod_ int egr al e 
use mod_creneau 
implicit none 

real : : yl , y2 

real , intrinsic exp 

yl=int egrale ( exp , . ,1.) 
y2=int egrale (creneau , — 1. ,1.) 
print * ,yl , y2 

end program prog_integrale 



module mod_creneau 




implicit none 




contains 




real function creneau( 


*) 


implicit none 




real , intent ( in ) 


X 


if (x<0.0.and.x>1.0) 


then 


creneau=0. 




else 




creneau = l. 




end i f 




end function creneau 




end module mod_creneau 
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5.4 Interface generique 

idee : on veut appliquer une meme procedure a des 
arguments de type differents. Exemple : abs(x) ren- 
voie \x\ si x est reel et V& 2 + b 2 si x = a + ib est 
complexe. 

- principe : dans un module, on ecrit plusieurs procedures 
adaptees a chaque type d' arguments, et on les re- 
groupe sous un meme nom "generique" . C'est le com- 
pilateur qui decide quelle procedure utiliser selon le 
type des arguments. 

exemple : on veut une fonction exponent ielle qui 
renvoie exp(x) si x est reel, et exp(z) = exp(a)(cos 6+ 
i sin b) si z = a + ib est complexe 
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module mod_exponent ielle 
implicit none 
interface exponent ielle 

module procedure rexp , zexp 
end interface 
contains 

function rexp(x) result(y) 
implicit none 
real , intent ( in ) x 
real y 

y=exp(x) 
end function rexp 

function zexp(z) result(y) 
implicit none 
complex, intent(in) z 
complex y 
real a , b 

a=real(z) ; b=aimag(z) 
y=cmplx (cos(b) , sin(b))*exp(a) 
end function zexp 

end module mod_exponent ielle 

methode : deux fonctions rexp et zexp stockees dans 
le module mod_ exponent ielle. La fonction generique 
exponent ielle est creee avec le bloc interface et 
l'instruction module procedure situee au tout debut 
du module. 
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- ensuite, dans un programme principal, on utilise la 
meme fonction exponent ielle, que l'argument soit 
reel ou complexe : 

program generique 

use mod_exponent ielle 
implicit none 

real x 
complex z 

x = 0. ; z=cmplx (0. ,3.1416) 

print*, 'exp(x)=' , exponent ielle ( x ) 
print*, 'exp(z)=' , exponent ielle ( z ) 



end program generique 
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5.5 Creation de nouveaux operateurs 

- exemple : on a une fonction rptens qui avec deux 
vecteurs reels v\ et en entree renvoie le produit 
tensoriel w — v\® 

w=r ptens ( vl , v2 ) 

- on voudrait utiliser cette fonction sous forme d'un 
operateur ptens : 

w=vl . ptens . v2 
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- on va utiliser un module, un bloc interface operator 



et l'instruction module procedure : 



module mod_ptens 




implicit none 




interface operator (. pt ens . ) 




module procedure rptens 




end interface 




contains 




function rpt ens ( vl , v2 ) result(w 


) 


implicit none 




real, dimension (:), intent ( in 


) : : vl , v2 


real , dimension ( size (vl ), size 


( v2 ) ) : : w 


integer :: i,j 




do i = l , size ( vl ) 




do i=l.size(vll 




w(i, j)=vl(i)*v2(j) 




end do 




end do 




end function rptens 




end module mod_ptens 





- les deux points autour de .ptens. sont obligatoires 
pour definir un operateur 

- dans le programme principal utilisant le module 
mod_ptens, on pourra utiliser l'operateur ptens ainsi : 



w=vl . ptens . v2 
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- on peut associer plusieurs fonctions a un meme ope- 
rateur (comme avec les interfaces generiques, sec- 
tion 5.4) : pour que .ptens. marche aussi avec les 
complexes, raj outer une fonction zptens dans le mo- 
dule, et ajouter le nom zptens a la suite de rptens 
dans l'instruction module procedure 

on peut aussi vouloir redefinir un operateur existant 
comme *,+,—,/ etc. Cela est possible, mais un peu 
plus complique (il faut que le nouvel operateur et 
l'ancien different par le type de leurs arguments). 
Dans l'exemple precedent, on ne pourrait pas redefinir 
* par l'operateur ptens car * est deja defini pour les 
tableaux. 
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- exemple de surcharge de + : pour concatener des 
chaines comme l'operateur //, en supprimant des 
blancs. Remarquer que + garde sa signification pour 
les autres types. 



program surcharge 

use mod_add_chaine 
implicit none 

character ( len =20) cl,c2,c 

cl='bon' ; c2='jour' 

c=cl+c2 

p r i n t * , c 

print *,' hello da'+' rling ' 
print* ,5 + 3 

end program surcharge 
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module mod_add_chaine 

implicit none 
interface operator(+) 

module procedure add_chaine 
end interface 

contains 

function add_chaine ( cl , c2 ) result(c) 
implicit none 

cha racter ( len =*) , intent(in) cl , c2 
cha racter ( len=len ( cl)+ len ( c2 ) ) c 

c=trim ( adj ustl (cl)) / / trim( adj ustl ( c2 ) ) 

end function add_chaine 

end module mod add chaine 
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5.6 Ressources privees, publiques, semi- 
privees 

- quand un programme utilise un module, il a acces 
a toutes les variables du module, et tous ses sous- 
programmes. Cela peut parfois provoquer des conflits 
avec les propres variables du programme si les noms 
sont identiques 

- exemple : un programmeur veut utiliser un module 
servant a calculer des normes de matrices avec la 
fonction norme pour son programme. Ce module con- 
tient beaucoup d'autres fonctions auxiliaires utilisees 
par la fonction norme. Le programmeur ne connait 
pas ces autres fonctions (il n'a acces qu'au fichier 
objet par exemple). Dans ce cas, les fonctions auxi- 
liaires ne devraient pas etre visibles a l'exterieur du 
module, car sinon, l'utilisateur n'aurait par exemple 
pas le droit de creer des fonctions portant le meme 
nom, alors qu'il ignore leur existence ! 

on peut alors utiliser les fonctions private (prive), 
public et les attributs correspondants private et public 

pour limiter 1' acces a certaines variables depuis l'ex- 
terieur du module 

- voir la reference [3] pour plus de details. 
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5.7 Point eurs 

- la structure de pointeur est fondamentale en langage 
C et C++. En fortran 90, une structure similaire 
existe, mais n'est pas d'un interet evident pour le 
calcul scientifique 

- dans la norme fortran 95, elle etait indispensable 
pour definir des tableaux allouables a l'interieur de 
types derives (comme par exemple pour definir des 
tableaux de tableaux). Mais depuis la norme fortran 
2000, cela n'est plus necessaire. 

- la notion de pointeur en fortran 90 ne sera done pas 
etudiee plus en detail dans ce cours 
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Liste de quelques fonctions intrinseques usuelles 

en Fortran 90 



adapte des cours de VIDRIS, disponibles sur 
http : //www. idris.Fr/data/cours/lang/f90/F90_cours .html 



On donne la signification de la commande, puis, dans un exemple, la syntaxe a utiliser, suivie 
apres le signe 'i— > ' de la valeur retournee en sortie. 



Fonctions intrinseques sur les entiers, reels, complexes, 
extensibles aux tableaux 

- ABS : retourne la valeur absolue de son argument. Pour un complexe, retourne sa norme : 

aA 2 + y 2 - 

ABS(-l) i-> 1 

ABS(-1.5) h-> 1.5; ABS((3.,4.)) ^ 5.0 

- ACOS : retourne l'arc cosinus en radians de son argument reel. 

AC0S(0. 54030231) h-> 1.0 

- AIMAG : retourne la partie imaginaire du complexe passe en argument. 

AIMAG((2. ,3.)) i-> 3.0 

- AINT : tronque le reel passe en argument 

AINT(2.783) i-> 2.0; AINT(-2.783) i-> -2.0 

- ANINT : retourne, sous forme d'un reel, l'entier le plus proche du reel transmis. 

ANINT(2.783) i-> 3.0 
ANINT(-2.783) i-> -3.0 

- ASIN : retourne Tare sinus en radians de son argument reel. 

ASIN(0. 84147098) i-> 1.0 

- ATAN : retourne l'arc tangente en radians de son argument reel. 

ATAN ( 1 .5574077) i-> 1.0 

- CEILING : retourne l'entier immediatement superieur au reel transmis en argument. 

CEILING (3.7) i-> 4 
CEILING(-3.7) i-> -3 

- CMPLX : retourne un complexe dont les parties reelle et imaginaire sont transmises en argument. 

CMPLX(-3.) i-> -3.0 + 0.i; CMPLX(2,4.) ^ 2.0 + 4.01 

- C0NJG : retourne le complexe conjugue de celui passe en argument. 

C0NJG(( -3.0 ,4.0)) ^ -3.0-4.01 

- COS : retourne le cosinus de Tangle passe en argument (exprime en radians). 

COS(l.O) i-> 0.54030231 
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- COSH : retourne le cosinus hyperbolique. 

COSH(l.O) i-> 1.5430806 

- DBLE : convertit en double precision l'argument transmis. 

- EXP : retourne l'exponentiel de l'argument transmis. 

EXP(l.O) i-> 2.7182818 

- FLOOR : retourne l'entier immediatement inferieur au reel transmis en argument. 

FL00R(3.7) i-> 3 ; FL00R(-3.7) i-> -4 

- INT : convertit en entier l'argument transmis. 

INT(-3.7) i-> -3 
INT(9. 1/4.0) i-> 2 

- LOG : retourne le logarithme neperien de l'argument transmis. 

L0G(2. 7182818) i-> 1.0 
LOG(IO.O) i-> 2.3025851 

- L0G10 : retourne le logarithme decimal de l'argument transmis. 

LOGIO(IO.O) i-> 1.0 
L0G10 (10.E10) i-> 11.0 

- MAX : retourne le maximum des nombres passes en argument. 

MAX( -9.0 ,7.0 ,2.0) i-> 7.0 

- MIN : retourne le minimum des nombres passes en argument. 

MIN( -9.0 ,7.0 ,2.0) h-> -9.0 

- MOD : retourne le reste de la division effectuee a l'aide des deux arguments fournis. 

MOD (3. 0,2.0) i-> 1.0 
M0D(-8,5) i-> -3 

- REAL : convertit en reel l'argument transmis. 

REAL (3) i-> 3.0 

- SIGN : retourne le nombre dont la valeur absolue est celle du premier argument et le signe celui 
du deuxieme. 

SIGN( -3.0,2.0) i-> 3.0 

- SIN : retourne le sinus de 1' angle passe en argument (exprime en radians). 

SIN(l.O) i-> 0.84147098 

- SINH : retourne le sinus hyperbolique. 

SINH(l.O) i-> 1.1752012 

- SQRT : retourne la racine carre de son argument. 

SQRT(5.0) ^ 2.236068010 

- TAN : retourne la tangente de Tangle passe en argument (exprime en radians). 

TAN(l.O) ^ 1.5574077 

- TANH : retourne la tangente hyperbolique. 

TANH(l.O) ^ 0.76159416 
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Fonctions intrinseques pour le temps de calcul 

CPU_TIME(t) (Norme 95) retourne dans le reel t le temps CPU en secondes (ou reel < 
si indisponible). Par difference entre deux appels, il permet d'evaluer la consommation CPU 
d'une section de code. 

call cpu_time ( t 1 ) 

call cpu_time(t2) 

print* , ' temps ^CPU^= ' ,t2-tl 

Fonctions intrinseques pour les chaines de caracteres 

AD JUSTL : cadre a gauche la chaine passee en argument : supprime les blancs en tete ; complete 
a droite par des blancs. 

ADJUSTL ( '_^Fortran ' ) ^ 'Fortran^' 

AD JUSTR : cadre a droite la chaine passee en argument : supprime les blancs en fin ; complete 
a gauche par des blancs. 

ADJUSTR ( ' Fortran^_ ' ) ^ ' _^Fortran ' 

LEN : retourne la longueur de la chaine de caracteres transmise en argument. 

CHARACTER ( len = 10) CH ; LEN(CH) ^ 10 

LEN_TRIM : retourne la longueur de la chaine de caracteres transmise en argument sans considerer 
les blancs de fin. 

LEN_TRIM( JDRTRAN^') i-> 9 
LEN_TRIM ( '_') i-> 

LGE : compare les deux chaines de caracteres transmises en argument : retourne .true, si la 
premiere chaine est superieure ou egale a la deuxieme, .false, sinon. 

LGE( 'MANET' , 'MONET') ^ .false. 

LGE ( 'MANET^Edouard ' , 'MANET' ) ^ . true . 

LGT : compare les deux chaines de caracteres transmises en argument : retourne .true, si la 
premiere chaine est superieure strictement a la deuxieme, .false, sinon. 

LGT ( 'MANET' , 'MANET' ) ^ .false. 

LLE : compare les deux chaines de caracteres transmises en argument : retourne .true, si la 
premiere chaine est inferieure ou egale a la deuxieme, .false, sinon. 

LLE ( 'MANET' , 'MONET') ^ .true. 
LLE ( 'MANET' , 'MANET' ) ^ . true . 

LLT : compare les deux chaines de caracteres transmises en argument : retourne .true, si la 
premiere chaine est inferieure strictement a la deuxieme, .false, sinon. 

LLT ( 'MANET' , 'MANET' ) ^ .false. 
TRIM : retourne la chaine de caracteres transmise debarrassee de ses blancs de fin. 

TRIM( 'PICASSO.uuuu') i-> 'PICASSO ' 
// : operateur pour concatener deux chaines. 

'bon' // 'jour' i— > 'bonjour' 
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4 Fonctions et subroutines intrinseques tableaux 



4.1 interrogation 

- SIZE (array [,dim] ) retourne la taille (ou l'etendue de la dimension indiquee via dim) du 
tableau passe en argument. 

integer, dimension (— 2 : 27 ,0 :49) :: t 
SIZE(t) ^ 1500 
SIZE(t , dim = l) ^ 30 

- MAXLOC(tab) et MINLOC(tab) retournent, pour le tableau tab, l'emplacement de l'element 
maximum/minimum (si le tableau est multidimensionel, le resultat est un vecteur dont les 
elements sont les indices de l'element dans tab) 

integer, dimension (— 2 : 27 ,0 :49) :: t 

integer, dimension ( 1 : 2 ) :: ind_elem_max , ind_elem_min 
ind_elem_max=maxloc ( t ) 
ind_elem_min=minloc ( t ) 

4.2 Reduction 

Selon que DIM est absent ou present, toutes ces fonctions retournent so it un scalaire soit un 
tableau de rang n-1 en designant par n le rang du tableau passe en premier argument. 

La plupart de ces expressions acceptent un "masque" en parametre, qui permet de n'appliquer 
les fonctions qu'a certains elements (voir le cours de 1'IDRIS de Fortran 95) 

- COUNT (verif) : verif est une expression logique sur deux tableaux conformes, count compte 
le nombre de fois ou cette expression est verifiee par deux elements 

| 1 3 5 | 

Soit A 

2 4 6 | 

Decompte global des valeurs de A > 2 : 

C0UNT(A>3) i-> 4 

- MAXVAL (array) et MINVAL (array) 

| 1 3 5 | 
Soit A | 

2 4 6 
maxval(A) i— > 6 
minval(A) i— > 1 

- PRODUCT (tab) renvoie le produit de tous les elements du tableau tab. 

PR0DUCT((/ 2,5,-6 /)) i-> -60 

- SUM (tab) renvoie la somme de tous les elements du tableau tab. 

SUM((/ 2,5,-6 /)) h- 1 

4.3 multiplication 

- D0T_PR0DUCT(vector_a, vector_b) retourne le produit scalaire des deux vecteurs passes en 
argument (produit hermitien si vector_a est complexe) 
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vl =(/ 2,-3,-1 /) ; v2 =(/ 6,3,3 /) 
DOT_PRODUCT(vl , v2) i-> 

MATMUL(matrix_a,matrix_b) effectue le produit matriciel de deux matrices ou d'une matrice 
et d'un vecteur passes en argument. Soit la matrice A = ^ !^ 3^ / ^ et le vecteur v — ^-4 j 

MATMUL ( A , V) 



4.4 transformation 

- TRANSPOSE (mat) renvoie la transposed de la matrice mat (tableau bidimensionnel) 
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